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^ I Abstract. This paper is concerned with the ergodic subspaces of the state spaces 

of isolated quantum systems. We prove a new ergodic theorem for closed quantum 
C*~) , systems which shows that the equilibrium state of the system takes the form of a 

grand canonical density matrix involving a complete commuting set of observables 
■ including the Hamiltonian. The result obtained, which is derived for a generic finite- 

^ ' dimensional quantum system, shows that the equilibrium state arising from unitary 

evolution is always expressible in the canonical form, without the consideration of a 
. system-bath decomposition. 
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PACS numbers: 05.30.-d, 05.30.Ch, 45.20.Jj 
Given the Hamiltonian H and the initial state |^o) of an isolated quantum system, 
what is the dynamic average 

((O)) = lim \ [\tp s \d\^ s )ds (1) 

>< 

of an observable O when the state \ipt) = e '^l^o) °f the system evolves unitarily? 
Is there an equilibrium density matrix p, with a thermodynamic characterisation, such 
that the average is given by ((O)) = tr(pO) ? 

In the case of a classical system, if the Hamiltonian evolution is ergodic, then the 
theorem of Koopman, von Neumann, and Birkhoff shows that the dynamic average can 
be replaced by a statistical average over a subspace of the phase space determined by the 
relevant conservation laws pQ. If the system consists of a large number of interacting 
particles, then the dynamic average is intractable, whereas the statistical average in 
many cases can be calculated. 

In the case of quantum systems, while the equilibrium properties of small 
subsystems of large systems have been studied extensively [2J, HI [5j El d El El HOI ITT] , 
less attention has been paid to the equilibrium states arising as a consequence of 
the unitary evolution of closed systems. The purpose of this paper is to investigate 
such systems and to derive rigorous results concerning (a) the dynamic averages of 
observables, and (b) the associated equilibrium states. 
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We consider an isolated quantum system based on a Hilbert space of dimension n+1, 
with a generic, nondegenerate Hamiltonian H (the degenerate case will be considered 
later). We write {-Ej}j=o,i,...,n for the energy eigenvalues, and = E-i — Ej for the 
eigenvalue differences. The normalised energy eigenstates will be denoted {|-Ei)}i=o,i,->ra> 
with the associated projection operators {ITj}i=o,i,...,n- We write \ipo) for the initial 
state, and {|V't)}o<i<oo for its unitary evolution under the influence of H. With these 
definitions at hand, the main result can be expressed as follows: 

Quantum ergodic theorem. The dynamic average of an observable O is given by 
{{O)) = tr(pO), where 

1 n 

and Z(j3, {/ij}) = tr exp(— j3H — X]™=2 A*^)- Here H together with {^} i=2 )Tl constitute 
a complete set of commuting observables. The effective inverse temperature (3 and 
chemical potentials {pi)i=2,..., n ar e given by the relations 
n 9S , dS 

P= dE> and ^ = W (3) 
where E = tr(pH), and Fj = tr(pFj). The entropy S = — tr(plnp) is given by 

n 

S = -J^ Pi In pi, (4) 

with pi = K^ol^)! 2 . 

The appearance of the grand canonical density matrix (j2J) is surprising, since this 
structure normally arises with the consideration of the equilibrium thermodynamics 
of a small system immersed in a thermal bath. Indeed, the canonical form p = 
exp(—j3H)/Z(f3) is known to appear in the case of a system in a thermal bath for 
an overwhelming majority of wave functions of the total system [TU1 [TT] . Equation ([2j) 
is a stronger result, valid in the case of a closed system, involving no approximations 
and no invocation of the thermodynamic limit. 

To determine the equilibrium states of a closed quantum system we need to identify 
the subspaces of the quantum state space over which a generic time evolution will 
exhibit ergodicity. The idea is that in general there are n conserved quantities arising in 
connection with unitary evolution in a Hilbert space of dimension These are given 

by the expectation values of n linearly independent observables that commute with the 
Hamiltonian, one of these being the Hamiltonian itself. Writing E for the expectation 
of H, we can then write {i^}i=2,...,n f° r the expectation values of the other members of 
the commuting set, which we denote by {i 7 i}j=2,...,n- By fixing the expectation values of 
these conserved quantities we are left with a set of n relative-phase degrees of freedom 
for the state vector that span the ergodic subspace of the state space associated with 
the given initial state. 

We shall show that the equilibrium state corresponds to a uniform distribution over 
the toroidal subspace of the quantum state space spanned by the relative phases. The 



Unitarity, ergodicity, and quantum thermodynamics 



3 



equilibrium distribution is characterised, in particular, by a density-of-states function 
Q, which acts as a measure of the size of the toroidal subspace. The associated density 
matrix p is given by the von Neumann-Luders state; that is to say, 

n 

P = ^2pi^i, (5) 

i=0 

where pj = \{^o\Ei)\ 2 . This might be surprising, since such a state arises most naturally 
in the context of measurement theory, where it describes the state of a system after an 
energy measurement has been performed. The result is consistent with the fact that the 
time average of the dynamics of the density matrix under unitary evolution is given by 
the von Neumann-Luders state. It follows that the dynamic average ([1]) of an arbitrary 
observable O is given by tr(pO). 

To identify the ergodic subspaces of the quantum state space, we first consider 
the example of a two- level system, with n = 1. The one-parameter family of states 
generated by unitary evolution can be written in the form 

|Vt> = cos^| J E 1 )+sin^e i ^ + ^|E ), (6) 

where < 9 < it and < < 27r. The pure state space has the geometry of a 
sphere, and unitary evolution gives rise to a rigid rotation of the sphere around the axis 
determined by the two energy eigenstates. Given the initial state the dynamical 
trajectory is the latitudinal circle on which |^o) lies. The circle is fixed by setting the 
initial energy E of the system, which is the only conserved quantity. Every point on the 
latitudinal circle is traversed by the dynamical trajectory, which makes this circle the 
ergodic subspace of the state space. The dynamic average of an observable can thus be 
replaced by the ensemble average with respect to a uniform distribution over the circle. 

To calculate the associated density of states we compute the weighted volume in 
the pure state manifold occupied by the states having the given property. In general, 
if we have a set of conserved quantities {Gj}j=i,..., m given by Gj = (ip t \G j\ipt) then the 
associated density of states is 

/m 
n sm&iW) -G}Wi» (7) 

3=1 

where the integration is over the space of pure states and dV^ is the associated volume 
element. The corresponding construction for classical systems is considered in [12J, 
where fl({Gj}) is referred to as a "substructure function". In the case of a two- 
level system the ergodic circle is chosen by fixing the expectation of the Hamiltonian: 
E = (H). In terms of the spherical coordinates (8, 0) of (EI), the constraint can be written 
intheform^-^cos 2 ^ = E-E . We thus integrate 5(cos 2 ~ 9-(E-E )/(E 1 -E )) 
over the pure state manifold. Since the volume element is dV^ = | sin 9d9d(f), the 
resulting density of states is 

n(E) = t {Eo<E<El} -^- r , (8) 

where ^-{a} denotes the indicator function: lm = 1 if A is true and Ij^j = otherwise. 



Unitarity, ergodicity, and quantum thermodynamics 



4 



We proceed to calculate the density of states for n — 2. In this case there are two 
conserved quantities: E = (H) and F = (F), where the observable F commutes with 
H, but is not of the form aH + bl. The calculation simplifies if we use an equivalent 
alternative set of constraints obtained by fixing the expectation values of two of the 
energy projectors, say po = (flo) and pi = (fli). It follows from the resolution of 
identity that p 2 = (n 2 ) = 1 — Po — Pi- The unitary trajectory can be written in the form 

\ip t ) = sin | $i cos 2 2 \E 2 ) + sin = 9 1 sin 

I^e^i+^l^^ + cosl^e^+^l^o), (9) 

and the two constants of motion are given by p = cos 2 | Qi and p 1 = sin 2 \ 9\ sin 2 \ 9 2 , 
which fix the variables 9i,9 2 . Therefore, under a generic unitary evolution the ergodic 
subspace of the quantum state space is the two-torus T 2 spanned by (f>i,(f) 2 . The 
density of states is obtained by integrating <5(cos 2 \ 9\ —p )S(sm 2 | 9i sin 2 | 9 2 — pi) over 
the pure state manifold, with the appropriate volume element, which in this case is 
dV = ^sinfl^l — cos^i) sin^ 2 d^id^ 2 d0id0 2 . Performing the relevant integration we 
find that Q(p ,Pi) = in the triangular region {0 < Po,Pi < 1} H {0 < p + p 1 < 1}, 
and vanishes otherwise. 

In the case of a general (n + l)-level system there are n conserved quantities 
associated with unitary dynamics. It follows that under a generic time evolution for 
which the eigenvalue differences {ouij} are incommensurate the typical ergodic subspace 
of the quantum state space is given by an n-torus T n . To calculate the density of states 
Q(po, ■ ■ ■ ,p n -i) we fix the constraints (fli) = pi for i — 0, . . . , n — 1, express these in 
terms of the coordinates (9i,(pi), and perform the constrained volume integral over the 
pure state manifold by using the volume element 

n 

dV = T n Y[ cos \ 9i sin 2 *" 1 \ 9^9^. (10) 
i=i 

The result is 

ft(P0,-",Pn-l)=7T n (11) 

in the hyper-triangular region {0 < po, ■ ■ ■ ,p n -i < 1} H {0 < po + • • ■ + p n -i < 1}) 
and il(po) • • • ,Pn-i) — otherwise. We see that irrespective of the Hilbert space 
dimensionality the density of states is constant in the hyper-triangular region, and is 
independent of the energy E and the conserved quantities {i 7 i}j= 2v ..,n- 

The analysis above leads to the following observation. Since for each n we have 
identified the ergodic subspaces of the state space, we are able to apply Birkhoff's 
theorem to conclude that the dynamic average of an observable can be replaced by the 
statistical average of the observable with respect to an equilibrium state given by a 
uniform distribution over the toroidal subspace T n . 

To compute the expectation of an observable O we determine the density matrix 
associated with the equilibrium distribution over the state space. We remark in this 
connection that the density matrix associated with a probability density function on 
the pure state manifold is the expectation of the pure-state projection operator with 
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respect to that density function. Now in the energy basis a pure-state projector can be 
expressed in the form 



Thus, the diagonal elements {pi} of the pure-state projector are real, whereas the off- 
diagonal elements contain phase factors. The equilibrium distribution has fixed values 
for the {pi} and a uniform distribution over the phase variables. It follows that if we 
take the average of the projector \ifj)(ip\ over the phases, the off-diagonal elements drop 
out and we are left with the von Neumann-Luders state (jSJ). 

The appearance of the von Neumann-Luders density matrix as the equilibrium state 
is consistent with the fact that the dynamic average of the density matrix is itself given 
by the von Neumann-Luders state. This can be seen as follows: 



In particular, we see that the timescale involved for the averaging to become effective 
is determined by the energy differences. We thus conclude that the dynamic average of 



density matrix, however, does not make the thermodynamic properties of the equilibrium 
state immediately apparent. We shall demonstrate, however, that in association with 
the conserved quantities (E, {Fi}) there is a corresponding system of conjugate variables 
(/3, {pi}) that can be given a consistent thermodynamic interpretation. In the case of 
the energy the conjugate variable has the interpretation of the inverse temperature. For 
the other observables the associated conjugate variables can be interpreted as chemical 
potentials. This suggests that the equilibrium state arising from unitarity and ergodicity 
might be of a grand canonical type. The conjugate variables are defined as follows. 
Writing (jSJ) for the density matrix associated with the toroidal subspace characterised 
by the conserved quantities (E, {Fi}) we have trp = 1, tv(pH) = E, and tr(pFk) = Fk- 
Let us define a family of n + 1 operators {G}j=o,i,...,n by setting G = 1, Gi = H, and 
{(jj}i=2,...,n = {Fi}i=2,..., n , writing {G}i=ox-,n for the corresponding expectation values 
with respect to p, so G — 1, G\ — E, and {Gi} i= 2,...,n — {^i}j=2,...,n- m other words, 
tr(pGi) — Gi for i = 0, 1, . . . , n. It follows from (jnj) that 




(12) 





(13) 



an observable O is given by tr(/)0), where p is given by This representation of the 




(14) 



i=0 
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Thus, writing = tr(fliGj) and defining hjk by Y!<j=o dv^jk = ^ifej we see that 

Pk = ^2Gjhjk, (15) 
i=o 

and therefore that 

n 

P = G i h jk^k- (16) 

i,fe=0 

To verify that hjk exists, we observe that if it did not, then there would exist a nonzero 
vector such that Y^j=o9ij^j = Oj but that would imply tr(fli Yll=o = f° r 

all z, and hence Yll=o A/'O = 0, contrary to the assumption that the Gj are linearly 
independent. 

Formula (fT6l) gives p as a function of E and {i^}- Therefore, writing S = — tr(pln p) 
for the entropy, we obtain an expression for S as a function of E and {i^}. The associated 
conjugate variables are then defined by the thermodynamic relation 

n 

dS = (3dE + J2^dE k , (17) 

fc=2 

where (3 is the effective inverse temperature and {pk} are the effective chemical 
potentials. This shows, on account of the linear independence of the observables, the 
equivalence of the specification of either (i) the initial state \i[)q) up to relative phases, 
(ii) the probabilities Pi = \{^ \Ei)\ 2 , (iii) the expectation values E and {-Fj}, or (iv) 
the conjugate variables (3 and {pi}. We can therefore investigate how the equilibrium 
density matrix ([5]) can be expressed either in terms of the extensive variables E and 
{Fi}, or in terms of the conjugate variables (3 and {pi}. 

For the various representations of the density matrix we consider first the example of 
the two-level system. In this case we solve the relations Po+pi = 1 and PqEq+piEi = E 
for the diagonal elements po,Pi of p, and obtain 

/ E 1 -E o \ 

m = El n E ° ■ (18) 

\ U Ei-Bo / 

Computing the entropy and using the relation dS = (3dE we can express the inverse 
temperature as a function of E. The result is 

By inverting this relation, we then obtain 

E W = ,-,» + *k • (20) 
Expression fl20|) is, however, the expectation of the energy with respect to the canonical 
density matrix. That is to say, (|18|) can be expressed in the form 

1 ( e"^ 
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where Z(/3) = e _/3So + e~ /3El . The important point here is that the inverse temperature 
(3 is not specified exogenously via the introduction of a heat bath. Rather, it is defined 
endogenously, through the specification of the energy of the equilibrium state associated 
with the given initial state. 

Let us now turn to the proof of the quantum ergodic theorem in the general case. 
It follows from ffl6l) that the entropy is given by 

n / n \ / n \ 

fc=0 \i=0 / \j=0 J 

Thus, defining 7$ = dS/ dGi by use of this expression, we find that 

n 

ii = -^2 hik 



k=0 

by (TT5]) . and hence 



j=Q 



^/i ife (lnp fe + l), (23) 



fc=0 



n n 

lnpi + 1 = -J2 9ijlj = ~ Yl tr A) 7j 
3=0 j=o 

n 

= - 7o - 7l tr (fL#) - T/tr (fl,^) . (24) 

Setting 71 = /? and {7i}j=2,...,n = {/^i}i=2,...,n, these relations are then sufficient to 
determine the diagonal elements {pi}i=o,...,n of the equilibrium density matrix in terms 
of the intensive variables, and we are led to the grand canonical ensemble (j2|) with 
the identification 70 = \nZ — 1. The effective inverse temperature, however, is not 
associated with an external heat bath, but rather is intrinsic to the system, and a similar 
remark applies to the effective chemical potentials. The fact that the conjugate variables 
are determined endogenously shows that our result does not require an assumption of 
entropy maximisation. 

In the case of a degenerate Hamiltonian, the ergodic subspace of the state space is 
contracted to a smaller torus T m C T n , where m + 1 is the number of distinct energy 
eigenvalues. This follows from the fact that since some of the eigenvalue differences LOij 
vanish, only m of the n relative phases for the unitary trajectory \ip t ) vary in time. As 
a consequence, we need only to consider m — 1 independent observables {i^} whose 
eigenspaces coincide with that of the Hamiltonian. In other words, there are only m 
terms, given by H and {i 7 i}j=2,...,m, in the exponent of (j2J) for the grand canonical density 
matrix. As an example consider the case of a three-dimensional Hilbert space where 
the energy eigenvalues are given by Eq, E±, and E\. The elements of the density matrix 
are p = (E x — E)/(E 1 — E ) and pi = p 2 = (E — E )/2(E 1 — E ). A short calculation 
making use of the relation dS* = /3dE then shows that 

m = E °l- m , + + \% . (25) 

which is evidently the expectation of H with respect to the canonical density matrix 
p = exp(-/3#)/trexp(-/3#). 
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A challenging open issue is to understand the implications of the quantum ergodic 
theorem for macroscopic systems. In the case of a large quantum system the energy 
spectrum of a typical model Hamiltonian is highly degenerate. As a consequence, 
the number of independent macro-observables H and {Fi} required for the exact 
specification of the equilibrium density matrix is significantly reduced. For real 
systems, however, due to the complexity of internal interactions one would expect the 
degeneracies in model Hamiltonians to split into closely located but distinct levels. 
Therefore, the specification of a small number of macro-variables will only provide an 
approximate description of the equilibrium state for real systems. On the other hand, 
if there are large clusters of observables with the property that in the equilibrium state 
defined by ([2]) the chemical potentials are approximately equal, then the resulting state 
can be adequately characterised by a small number of macro-variables, and thus can 
be regarded as effectively classical. It is interesting in this connection to contrast the 
results obtained here for quantum systems with the corresponding results for strictly 
classical systems: While ergodicity is generic for quantum systems, classically it is 
exceptional [13J . The fact that the characterisation of the equilibrium state of a quantum 
system is simpler, and that the equilibrium distribution can be derived dynamically by 
use of an ergodicity argument, might be related to the special structures of energy 
surfaces in quantum phase spaces 

DCB acknowledges support from The Royal Society. 
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